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We study the potential of the TESLA hnear cohider operated at a 
center-of-mass energy of 500 to 1000 GeV for the measurement of the 
neutral Higgs boson properties within the framework of the MSSM. 

r~| , The process of associated Higgs boson production with subsequent 

decays of Higgs bosons into b-quark and r-lepton pairs is considered. 
An integrated luminosity of 500 fb~^ is assumed at each energy. The 
Higgs boson masses and production cross sections are measured by 

5^ \ reconstructing the bbbb and bbr+T~ final states. The precision of 

these measurements is evaluated in dependence of the Higgs boson 
masses. Under the assumed experiinental conditions a statistical ac- 
curacy ranging from 0.1 to 1.0 GeV is achievable on the Higgs boson 
mass. The topological cross section (j{e^e~ — > HA -^ bbbb) can be 
determined with the relative precision of 1.5 - 6.6 % and cross sec- 
tions a{e&e~ — > HA — > bbr^r") and cr(e"'"e~ — > HA -^ r"'"r~bb) with 
precision of 4 - 30 %. Constraints on the Higgs boson widths can be 
set exploiting bbr+r" channel. The ba discovery limit corresponds 
to the Higgs mass of around 385 GeV for the degenerate Higgs boson 
masses in the HA -^ bbbb channel at y/s = 800 GeV with integrated 
luminosity of 500 fb~^. The potential of the Higgs mass determination 
for the benchmark point SPS la for the process HA — > bbbb at y/s = 
1 TeV and luminosity 1000 fb~^ is investigated. 



1 Introduction 

Elucidating the mechanism of electroweak symmetry breaking will be one of 
the central tasks at a future linear e"'"e~ collider. In the minimal Standard 
Model (SM) breaking of electroweak symmetry is provided by one Higgs 
doublet, which introduces an additional spin-0 particle, the Higgs boson [1]. 
Among models with an extended Higgs sector, two Higgs doublet models 
(2HDM) [2] are of particular interest, since this structure of the Higgs sector 
is required in the Minimal Supersymmetric Standard Model (MSSM) [3], 
the most seriously considered extension of the SM. The Higgs sector of the 
CP-conserving MSSM comprises the following physical states: two CP-even 
Higgs bosons, the lighter of which is denoted h and the heavier H, one CP-odd 
Higgs boson A and two charged bosons H^. Two crucial parameters, which 
influence the Higgs boson phenomenology, are tan/3, the ratio of vacuum 
expectation values of the two Higgs doublets, and the angle a, which controls 
mixing in the CP-even sector. 

The cross sections of the Higgs-strahlung and fusion processes, involving 
CP-even Higgs bosons, 

e+e" -^ hZ(HZ), e+e" -^ h(H)z/ei^e, e+e" -^ h(H)e+e-, (1) 

are defined in terms of the strength of the Higgs couplings to weak bosons 
at tree level. These depend on the parameters a and /3 in the following way: 

ghzz.hww ~ sin(/5 - a), gHzz,Hww ~ cos(/5 - a). (2) 

In the MSSM, the set of the tree level couplings between neutral Higgs par- 
ticles and weak bosons is extended by two additional couplings, 

ghAz ~ cos(/5-a), gHAz ~ sin(/5-a). (3) 

As a consequence, in e"'"e~ collisions, the Higgs-strahlung and fusion pro- 
cesses will be complemented by the associated Higgs boson pair production 
mechanisms: 

e+e~ -^ hA(HA). (4) 

To fully explore the MSSM Higgs sector, the processes (4) must be stud- 
ied, since they directly probe ghAz and gHAz trilinear couplings, provide access 
to the A boson properties and supplement the study of the CP-even Higgs 
boson properties via Higgs-strahlung and fusion mechanisms. 

The study of the Higgs pair production process is motivated by a so- 
called decoupling limit of the MSSM in which the h boson approaches the 
properties of the SM Higgs boson. The closer MSSM scenario moves towards 
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the decoupling limit the more difficult it becomes to distinguish the Higgs 
sector from the SM. In such a scenario the detection of heavier neutral Higgs 
bosons would be crucial for establishing an extended Higgs sector. The de- 
coupling limit is approached relatively fast for large values of the H and A 
boson masses, ttia? "^h > 200 GeV, in a large portion of the MSSM parameter 
space. The distinct feature of this scenario is vanishing coupling of heavy 
CP-even Higgs boson to weak bosons, cos(/3 — a) -^ 0. As a consequence, 
the H boson production via the fusion and Higgs-strahlung processes is sig- 
nificantly suppressed, whereas the cross section of the e^e~ — >■ HA process 
reaches its maximal value making associated heavy Higgs pair production a 
promising channel for the detection of the H and A bosons at a future lin- 
ear e"''e~ collider. It should be also emphasized that in the decoupling limit 
the H and A bosons are almost degenerate in mass and have similar decay 
properties^. 

In this note the study of Higgs boson pair production at a future linear 
e"''e~ collider is presented. This analysis extends the previous studies [4, 5]. 
The analyzed topologies include the bbbb and bbr"'^r~ final states. In the 
following we will assume the e"'"e" -^ HA production process. However, this 
assumption does not restrict the generality of our study and the analyses 
developed in this note are applicable to the process e"'"e~ -^ hA as well as 
to the processes of associated Higgs pair production in CP-violating MSSM 
scenarios where the neutral Higgs bosons are no longer CP eigenstates [6]. 

This note is organized as follows. In Section 2 the assumed experimental 
conditions and main features of the detector are discussed. Section 3 de- 
scribes the simulation of the signal and the relevant background processes. 
In Section 4 the analysis tools, employed in our study, are briefly reviewed. 
Section 5 describes the analysis of the bbbb and bbr"*~r~ final states. In 
Section 6, we discuss the prospects of measuring the Higgs boson properties 
using the bbbb and bbr+r" topologies. In Appendix the interpretation of 
this analysis is done for the SPS la benchmark point. 

2 Experimental Conditions and Detector 
Simulations 

The study is performed for a linear collider operated at a center-of-mass 
energies of 500 GeV to 1000 GeV and an event sample corresponding to an 
integrated luminosity of 500 fb~^. 

The detector used in the simulation follows the proposal for the TESLA 
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collider presented in the TDR [4]. Studies were also performed in Asia and 
North America [7, 8]. The interaction region is surrounded by a tracker 
consisting of a multi-layered pixel micro-vertex detector (VTX) as the in- 
nermost part and a large time projection chamber (TPC) supplemented 
with additional silicon strip detector in the forward region (FTD) and two 
bands between VTX and TPC. In radial direction follow an electromagnetic 
calorimeter, a hadron calorimeter, the coils of a superconducting magnet 
and an instrumented iron flux return yoke. The solenoidal magnetic field is 
4 Tesla. 

The envisaged impact parameter resolution provided by the vertex detec- 
tor is 

SiIPr^,z) < 5/im © ^ . 3/, / , (5) 

psm' 6 

where p is the track momentum and 9 is the polar angle. The overall central 
tracker momentum resolution is 

ct(-) <5-10~^ [GeV/c]-\ (6) 

where pt is the transverse momentum. To enable efficient separation of 
neutral and charged particles, highly granular electromagnetic and hadronic 
calorimeters are foreseen for the TESLA detector. Their energy resolutions 
are: 

""^' - '°^° © 1%, ^ = ^©4%, (7) 
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where Eei and Eh are the energies measured for electrons and hadrons in the 
corresponding calorimeters. The polar angle coverage of the central tracker 
maintaining the resolution is |cos6'| < 0.85, above this range the tracking 
resolution deteriorates. The electromagnetic and hadron calorimeters cover 
I cos^l < 0.996 maintaining the resolution over the whole angular range. 

The detector response is simulated with parametric fast simulation pro- 
gram SIMDET [9]. Beamstrahlung is accounted for using CIRCE [10]. 

The event reconstruction is done in terms of energy flow objects. First, 
tracks are measured with tracking system and associated to calorimetric clus- 
ters to define charged energy flow objects (electrons, muons and charged 
hadrons). Since momentum measurement with tracking system is much 
more accurate than direction and energy measurements with calorimeters, 
the tracking information is used for estimation of the four-momentum of the 
charged objects. Calorimetric clusters with no associated tracks are regarded 
as neutral energy flow objects (photons and neutral hadrons). Measurements 



of the four-moinentum of neutral objects are solely based on the calorimetric 
information. 

3 Physics Processes and Monte Carlo 
Samples 

Samples of e"'"e~ -^ HA events are generated for several Higgs boson mass 
hypotheses with PYTHIA 6.2 [11], including initial state radiation. The tree 
level cross section of the HA production is related to the Higgs-strahlung 
cross section in the SM, af^, in the following way [2]: 

0"HA = gHAZ-^HAO-RZ ' (8) 

where 

,3/2 
^HA = ,,, 7^ (9) 

Al/|[12m|/s + AHz] 

accounts for the correct suppression of the P-wave cross section near the kine- 
matic threshold. The quantity Ajj = [1 — (mj+mj)^/s][l — (mj— mj)^/s] is the 
usual momentum factor of the two particle phase space. The Higgs-strahlung 
cross section in the SM is computed using the program HPROD [12]. 

For the case of maximal allowed gHAz coupling, sin(/3 — a) = 1, which 
is assumed for these studies, signal cross sections are given in Table 1. The 
branching fractions of the Higgs bosons into b-quarks and r-lepton pairs are 
set to their typical values in the MSSM : Br(H -^ bb) = Br(A -^ bb) = 90%, 
Br(H -^ tt) = Br(A -^ tt) = 10%. 

For the background estimation, the following processes are generated us- 
ing PYTHIA 6.2: e+e" -^ WW, e+e" ^ ZZ, e+e" -^ qq^, e+e" -^ Ze+e" 
and e"''e~ — > Weu. The cross sections for the most important background 
processes are given in Table 2. 



4 Analysis tools 

4.1 B-tagging 

Identification of b-quarks plays a crucial role in this analysis. Efficient tag- 
ging of jets containing heavy flavour hadrons will be achieved with a highly 
granular micro-vertex detector, allowing for precise reconstruction of track 
parameters in the vicinity of the primary interaction point. 
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The procedure of tagging b-jets exploits information from the single track 
as well as secondary vertex information. Secondary vertices are searched 
within jets using the package ZVTOP [13] developed for the SLD experiment, 
where tracks are described as probability tubes and seed vertices are defined 
as regions where these tubes overlap. Afterwards an attempt is made to 
assign additional tracks by an iterative procedure. For each found vertex the 
invariant total mass and momentum are calculated from the four-momenta of 
energy flow objects assigned to the vertex. Three dimensional decay length 
and decay length significance are also computed. 

For a jet flavour separation a neural network is developed. The detailed 
description of the neural network implementation for the full simulation of 
the TESLA detector with the BRAHMS program [14] can be found in [15]. 
Three different neural networks are introduced for jets with no, one and more 
than one secondary vertices found. For jets without secondary vertex, the 
impact parameter joint probability [16] and the two highest impact parameter 
significances are used as inputs for the neural network. If jet contains at least 
one secondary vertex, an additional information, including vertex masses, 
momenta, decay lengths and decay length significances, is fed into the neural 
network. 

The neural networks are trained on event samples simulated with SIMDET 
using the same variables and jet classification as in BRAHMS. For the anal- 
ysis, a jet-wise tag, referred hereafter as jet b-tag variable, is used. For a jet 
with neural network output x it is defined as 

where fb{x) and fudsc{x) are probability density functions of neural network 
output in a samples of b- and udsc-jets, respectively. Tagging of c-jets 
proceeds in a similar way. 

Fig. 1 shows b-tag and c-tag purity versus efficiency curves for Z ^ qq 
events simulated at a center-of-mass energy of ^/s = 91.2 GeV. Results ob- 
tained with SIMDET and BRAHMS are compared in this figure. The c-tag 
performance agrees within 5% over the entire range of efficiency. Some dis- 
crepancy in modelling of the b-tag is observed in the region of high efficiency 
and low purity. This discrepancy occurs due to not adequate modelling of 
the resolution tails in the impact parameter joint probability distribution by 
SIMDET. However, both bbbb and bbr^r" analyses impose b-tag require- 
ments strong enough not to be sensitive to this discrepancy. Some flavour 
tag related systematic studies have been performed. It has been found that 
the b-tag performance is nearly independent of the center-of-mass energy in 
the range from i/s = 91.2 GeV to -y/s = 500 GeV. Furthermore, possible 



changes in micro-vertex detector configuration are found to liave marginal 
impact on tlie b-tag performance. For example, removing innermost silicon 
layer changes the selection efficiency by not more than 5%. These studies 
confirm the stability of b-tagging. However, c-tag and e.g. b-quark charge 
tagging are depending on the innermost layer. 

4.2 Kinematic Fit 

The mass resolution of the reconstructed Higgs bosons is improved by means 
of a kinematic fit. In the bbbb analysis, conservation of four- momentum 
is required, leading in total to four constraints. In the case of the bbr''"r~ 
topology, the energies of r-leptons are not measurable due to undetectable 
neutrinos. However, the flight direction of highly boosted r-leptons can be 
approximated with good accuracy by the direction of the total momentum of 
its visible decay products. Exploiting this approximation together with the 
the four-momentum conservation requirement, we arrive at a two constraint 
fit. 

The jet angular and momentum resolution functions and the angular 
resolution of r-leptons are derived from Monte Carlo studies. The kinematic 
fit is performed using the code developed by DELPHI [17]. For the samples of 
the HA ^ bbbb and HA -^ bbr+r~ events with (mn, tua) = (150,100) GeV, 
the performance of the kinematic fit is illustrated in Fig. 2 and 3, respectively. 

5 Analysis Procedures 

5.1 The HA -^ bbbb Channel 

Events of the bbbb topology are characterized by four high multiplicity 
hadronic jets, containing decay products of b-hadrons. A cut-based tech- 
nique is employed to separate signal from background. Selection criteria are 
optimized separately for 500 GeV and 800 GeV center-of-mass energies. Each 
event is required to pass the following cuts: 

1. Total energy deposited in the detector, the visible energy Evis, must be 
greater than 340 GeV (600 GeV) for 500 GeV (800 GeV) center-of-mass 
energies. 

2. Each event is forced into four jets using the DURHAM algorithm [18] 
and the number of tracks per jet is required to be greater than three. 



3. To separate centrally produced H and A bosons from the WW and ZZ 
events, peaking in forward/backward direction, we apply a cut on the 
polar angle of the thrust vector [11], | cos^^-l < 0.8. 

4. Further suppression of the WW and ZZ backgrounds is achieved by 
requiring the event thrust value to be less than 0.88. 

5. Two fermion background is suppressed by applying a cut on the 
DURHAM jet resolution parameter, for which the event changes from 
four to three jets, \0giQys4 > —2.9. 

6. High multiplicity six-jet events originating from the e"''e~ -^ tt pro- 
duction are reduced by requiring the number of enflow (energy flow) 
objects in the event to be less than 130. This cut is applied only at ^/s 
= 500 GeV. 

7. The background from e"''e~ — > tt events is further reduced by applying 
a cut on the jet resolution parameter, for which the event changes from 
six to five jets, logioZ/se < -3.1 (-2.8) at ^i = 500 GeV (800 GeV). 

8. Finally, we make use of the b-tag information to enhance the purity 
of the selected event sample. First, the b-tag variable for each jet is 
calculated as described in Section 4. The four b-tag variables are sorted 
in descending order, Bi > B2 > B^ > B4. Two quantities B12, B^^ are 
then defined as 

7-) -D1-D2 

'^^ B,B2 + {l-B^){l-B2y 
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5354 + (l-53)(l-54) 

The value of -B12 has to be greater than 0.75 (0.6) at ^i = 500 GeV (800 
GeV). The value of 534 is required to be greater than 0.05 independent 
of the center-of-mass energy. 

The numbers of expected signal and background events, retained after selec- 
tion, and signal efficiency for the example points with {mi[,mx) = (150,100) 
GeV at v^ = 500 GeVand (mH,mA) = (300,250) at v^ = 800 GeVare pre- 
sented in Tables 3, 4. Fig. 4 and 5 show the distributions of the selection 
variables for the different sources of background and for the signal at a/s = 
500 GeV. The cutflow table for this case shows that such cuts like those on 
visible energy, number of tracks per jet, polar angle on the thrust vector and 



b-tag cuts are very effective to separate signal from the background. Espe- 
cially should be stressed that after all kinematic selection i.e. before b-tag 
cuts the signal efficiency is 54% while the background is drastically reduced 
(1.6% from the initial value). After b-tag cuts only 1.3% of the background 
events, left after kinematic selection, remains, while the signal efficiency is 
43%. This confirms, that the new tool for b-tagging is very powerful. The 
signal efficiencies, number of signal events and total background for different 
Higgs boson masses after selection cuts are given in Table 5. 

Events accepted in the final sample are subjected to a 4C kinematic fit. 
For each of the three possible di-jet pairings, the di-jet mass sum and the 
di-jet mass difference are reconstructed. For the example point the mass sum 
and the mass difference presented in Fig. 2 show that due to the kinematic 
fit the mass reconstruction resolution can be improved significantly. 

5.2 The HA -^ hhr+r- ,T+T-hh Channels 

The signature of the bbr'^r^ final state consists of two high-multiplicity 
hadronic jets enriched with the decay products of b-hadrons and two low- 
multiplicity jets initiated by r-leptons. 

Initially, a preselection of high-multiplicity events compatible with the 
four-jet topology is applied. Events are required to fulfill the following crite- 
ria. 

1. The leptonic and two photon events are rejected by applying lower cut 
on the number of energy fiow objects, Ncnflow > 30. 

2. Genuine six-jet events, resulting from the e"'"e~ -^ tt production and 
characterized by very large particle multiplicity, are partially suppressed 
by an upper cut on the number of energy fiow objects, Ncnflow < 120. 

3. The visible energy E^is must be greater than 250 GeV (400 GeV) at 
^/s = 500 GeV (800 GeV). This cut suppresses background processes 
characterized by large missing energy: ZZ -^ qqz/i^ events, radiative re- 
turns to the Z resonance with a photon escaping into the beam pipe 
and e"'"e~ — > Wez/, e^e~ — > Ze~'"e~ processes with electrons escaping un- 
detected in the forward direction. An upper cut Evis < 500 GeV (760 
GeV) at ^/s = 500 GeV (800 GeV) allows partial rejection of the fully 
hadronic WW -^ qqqq and ZZ — > qqqq events without significant loss 
in the signal efficiency. 

4. The event thrust value has to be less than 0.95. This cut disentangles 
more spheric HA events from the WW, ZZ and qq backgrounds in which 
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the high boost of the final state bosons/quarks results in large values 
of the event thrust. 

5. To separate centrally produced H and A bosons from the WW and ZZ 
processes in which weak bosons tend to peak in the forward/backward 
directions, a cut on the polar angle of the thrust vector, | cos 6*^1 < 0.9, 
is applied. 

6. The event is discarded if it contains an e+e~ or fi~fi~^ pair with invariant 
mass compatible with the Z mass within 5 GeV. This cut suppresses 
the ZZ — > qqi'^£~ background. 

Events satisfying these criteria are resolved into four-jet topology and r- 
lepton identification is performed. Two identification categories for r-leptons 
are introduced. A low multiplicity jet is assigned to the first category if it 
fulfills one of the following tight conditions: 

• A jet contains only one charged track, identified as an electron or muon. 
No neutral clusters are assigned to the jet. 

• A jet contains one charged track, identified as hadron, and not more 
than four neutral clusters (one-prong hadronic decay of r-lepton). 

• A jet contains three charged tracks with unit total charge and has 
not more than two neutral calorimetric clusters (three-prong hadronic 
decay of r-lepton) . 

To enhance the acceptance of signal events, the second identification category, 
imposing looser criteria, is introduced. According to these criteria, a jet is 
identified as r-lepton if it has not more than two charged tracks and less than 
five neutral clusters. 

An event is accepted if at least one of the low multiplicity jets satisfies 
the tight identification criteria. Signal efficiencies and numbers of signal and 
background events remaining after preselection are given in Table 6. 

In the next step, events are selected into the final sample by means of 
a binned likelihood technique. The signal likelihood, Lha is built from the 
following variables: 

• event thrust, T; 

• the polar angle of the thrust vector, cos 9t', 

• minimal opening angle between any two of the four jets; 



• b-tag variables of the two hadronic jets, Bi and B2 (hadronic jets are 
ordered by their energy); 

• missing energy; 

• number of energy flow objects, Nenflow 

Fig. 6 and 7 show the distributions of these variables. The distribution of 
Lha is shown in Fig. 8 for (mH,mA) = (150,100) GeV at ^i = 500 GeV. 
An event is accepted in the final sample if the value of Lha is greater than 
certain threshold optimized separately for each Higgs boson mass hypothesis 
considered to yield maximal significance of the signal. The number of signal 
and background events and the signal efficiencies after final selection are given 
in Table 7. Signal efficiency ranges between 30 and 45%, and the number of 
background events contributing to the final sample between 2000 and 6000, 
depending on the Higgs mass. 

For events accepted in the final sample, the di-jet and di-tau invariant 
masses are reconstructed exploiting the 2C kinematic fit, described in Section 
4. 

6 Results 

In the final step of the analysis, the spectra of the di-jet mass sum and 
difference, obtained in the HA — *■ bbbb channel, and distributions of the di- 
jet and di-tau masses, reconstructed in the HA -^ bbr'^r", r"''r~bb channels, 
are used to determine Higgs boson properties. 

6.1 Cross section and Mass 

First, the analysis is performed assuming that the natural widths of the 
Higgs bosons are small compared to the detector resolution. As an example 
Fig. 9 shows the distributions of the di-jet mass sum and di-jet mass difference 
obtained after selection cuts and kinematic fit in the HA — ^ bbbb channel for 
the Higgs boson mass hypothesis of (mH,mA) = (300,250) GeV at y/s = 800 
GeV. Three entries per event contribute to these distributions, corresponding 
to three possibilities to pair jets in the four-jet events. Two entries form 
a so-called combinatorial background. Fig. 10 demonstrates the final di- 
jet mass sum and di-jet mass difference after the cut on the di-jet mass 
difference sum, respectively, as indicated by arrows in Fig. 9. The signal 
efficiencies, number of signal events and total background for different Higgs 
boson masses after cuts on di-jet mass sum and difference are given in Table 5. 
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The mass distributions are fitted separately with the signal normalization, 
the sum and the difference of the Higgs boson masses as free parameters. The 
shapes of the signal distributions are parametrised using high statistics signal 
samples. Background is approximated with a polynomial function. From the 
fit, the sum and the difference of the Higgs boson masses and errors on these 
quantities are determined. The di-jet mass sum and difference are found to 
be weakly correlated quantities. Hence, the statistical errors on the masses 
of the Higgs bosons can be calculated as follows: 



where (5S and 6 A are statistical uncertainties in determination of the Higgs 
boson mass sum and difference. The error on the topological cross section is 
calculated as 

6a /a = ^/Ns + Nb/Ns. 

The notations Ns and A'^^ stand for the number of background and signal 
entries within the windows in the di-jet mass sum and di-jet mass difference 
distributions, where the signal is accumulated. The boundaries of these win- 
dows are optimized to yield a minimal relative error on the topological cross 
section. 

In the HA -^ bbr+r", r+r~bb channels, in the case of large mass splitting 
between Higgs bosons, the signal exhibits itself as two peaks in the recon- 
structed di-jet and di-tau mass spectra as illustrated in Fig. 11. The spectra 
are fitted simultaneously with the superposition of background and signal 
distributions. The fit is performed with four free parameters : two normaliza- 
tion factors for the HA — > bbr+r" and HA -^ r"'~r~bb samples and two Higgs 
boson masses. Background distribution is assumed to be well measured and 
therefore fixed. The shapes of signal peaks are parametrised using high statis- 
tics MC samples. When H and A are degenerate in mass, signal distributions 
overlap and the contributions from the HA — »■ bbr^r" and HA — *■ r"'"r~bb 
samples cannot be disentangled. To stabilize the fitting procedure we assume 
that the numbers of the HA -^ bbr"''r~ and HA -^ r^r^bb events, selected 
into the final sample, are equal. This assumption is validated by the proper- 
ties of the decoupling limit. As was mentioned in the introductory section, 
in the decoupling limit the H and A bosons have almost identical proper- 
ties: mn ~ m-A, Br(H -^ bb) ^ Br(A — > bb), Br(H -^ tt) ^ Br(A -^ tt), 
Th ~ Ta- Consequently, the HA -^ bbr+r" and HA -^ r+r^bb events are 
expected to contribute equally to the final selected sample. This reduces the 
number of free parameters to three : a common normalization factor for the 
two signal samples and two Higgs boson masses. As an example. Fig. 12 
shows the result of the fit for (mH,mA) = (300,300) GeV. 
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Table 8 summarizes the statistical accuracy of the mass measurements 
for the HA -^ bbbb and HA -^ bbr^r^, r+r^bb channels and their combi- 
nation. The statistical accuracy of the combined mass measurement is ~ 
100 MeV for the Higgs pair production far above kinematic threshold and 
degrades to ~ 1 GeV with aproaching kinematic limit. At the same time the 
statistical errors for the ^/s = 800 GeV are in general twice as bigger than 
those for the 500 GeV for the both channels. The combination of the two 
channels helps to improve the mass determination accuracy. 

The statistical uncertainty on the topological cross sections is reported 
in Table 9. In the HA — *> bbbb channel, the topological cross sections can be 
measured with relative precision between 1.5 and 6.6% while in the 
HA -^ bbr+r^, r+r^bb channels between 5 to 30%. 

A big part of SUSY parameter space leads to degenerate H and A Higgs 
boson masses. For this case the discovery significance as a function of tuh 
{=mA) (Fig. 13) is calculated for the HA -^ bbbb channel at ^/s = 800 GeV. 
Approaching the kinematic limit, the significance drops below 5a between 
380 GeV to 390. The whole range of the significances for the Higgs masses 
is from 28.2 to 3.4. The 5a discovery reach and 77^ as a function of rriH 
(=m^) where 77^ is assumed e"''e~ -^ HA cross section relative to that for 
sin^(/? — a) = 1 (as shown in Fig. 13). 

6.2 Width 

We utilize the HA — ^ bbr^r^, r+r^bb channels to measure the Higgs boson 
widths. Fig. 14 illustrates the impact of the natural widths of the Higgs 
bosons on the reconstructed di-jet and di-tau mass spectra. We consider 
separately two cases : 

1. large mass splitting between H and A, 

2. mass degeneracy. 

In the first case, the strategy of the measurement is the following. The masses 
of Higgs bosons are fixed to the values measured in the HA -^ bbbb chan- 
nel. The reconstructed di-jet and di-tau mass line shapes are parametrised 
as a function of Fh and Fa using MC samples of the HA — *> bbr"'"r~ and 
HA —* r'^r^bb events generated with different Higgs boson widths. A simul- 
taneous log-likelihood fit of the parametrised mass line shapes to the actual 
reconstructed di-jet and di-tau mass distributions is performed with four free 
parameters : Fh, Fa and two normalization factors for the HA -^ bbr"'"r~ 
and HA — > r"'"r~bb samples. Fig. 15 shows the dependence of the log- 
likelihood on the probed Higgs boson widths for the case of Fh = Fa = 
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5 GeV and {mn,mA) = (200,150) GeV at v^ = 500 GeV. The dependence is 
shown for the case when only di-jet mass or di-tau mass spectrum are fitted 
and for the simultaneous fit to the two distributions. The feasible accura- 
cies of the width measurements for three representative Higgs boson mass 
hypotheses are given in Table 10. The precision on the width determination 
ranges from 600 MeV to 4.0 GeV. 

In the case of very close Higgs boson masses, the broadening of the signal 
peak in the mass distributions can be caused by two factors : the finite 
natural widths of the Higgs bosons and non-zero mass difference, 6M = mu — 
rriA- The width measurement in this case proceeds as follows. The sum of the 
Higgs boson masses is fixed to the value measured in the HA — > bbbb channel, 
mn + '^A = 2M. From the HA -^ bbr+r^ and HA -^ r+r~bb MC samples 
with Fh = Fa = GeV, the experimental mass resolution functions Rbb^rrix) 
are derived, which approximate the reconstructed di-jet and di-tau mass 
distributions in the limit of infinitely small Higgs widths. To ensure stability 
of the fitting procedure, the H and A boson are assumed to have identical 
decay properties, Br(H -^ bb) = Br (A -^ bb), Br(H -^ tt) = Br (A — > tt), 
Th = Ta = F, as expected in the decoupling limit. The distributions of the 
invariant di-jet and di-tau masses are fitted simultaneously with the functions 

Fbb,rr{x) = Bbb,rr{x) + N ■ BW {x, M ± 5 M / 2,V) ® Rbb,rr{x) + 

A^ ■ BW{x, M T 5M/2, F) ® RbhM^) 

The fitting functions consist of three terms. The first one, Bbb^Trix)-, describes 
the background spectrum. The other two describe two signal distributions 
and represent the convolution of the Breit-Wigner function 

tMV 
BWix M T) = — 

^' '^ (X2 - M2)2 + M2F2 

with the experimental resolution functions Rbb^rrix)- The fit is performed 
with the common normalization factor, A^, the Higgs boson width, F, and 
the Higgs boson mass difference, 5M, as free parameters. The Higgs width 
errors obtained from the fitting procedure are presented in Table 11 for three 
representative Higgs mass hypotheses. The precision on the width determi- 
nation ranges from 1 GeV to 4.5 GeV. 

7 Conclusion 

We examined the potential of a future linear collider detector for the determi- 
nation of the MSSM Higgs boson properties exploiting associated Higgs boson 
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pair production followed by the Higgs decays to b-quarks and r-leptons. It is 
shown that combining the HA -^ bbbb and HA -^ bbr"''r~, r+r~bb channels 
the Higgs boson masses can be measured with an accuracy of up to several 
hundred MeV for Higgs pair production far above the kinematic threshold. 
The precision deteriorates to about 1 GeV with approaching the kinematic 
threshold. The topological cross section o"(e"'"e~ — ^ HA -^ bbbb) can be mea- 
sured with a relative precision varying between 1.5 and 6.6%. The error on 
the topological cross sections cr(e'^e~ — ^ HA -^ bbr"'"r~) and 
a{e'^e~ — > HA -^ r'^r~bb) is estimated to range between 4 and 30%. Mod- 
erate constraints on the Higgs boson widths can be set from the analysis of 
the reconstructed di-jet and di-tau mass spectra in the HA -^ bbr+r~ and 
HA -^ r+r^bb channels. 

The 5cr discovery limit corresponds to the Higgs mass of around 385 GeV 
for the degenerate Higgs boson masses in the HA -^ bbbb channel at ^/s = 
800 GeV with integrated luminosity of 500 fb~^. 
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A SPS la 

The present analysis is applied to one of the so-called benchmark points 
SPS la for SUSY searches [19]. SPS 1 is a typical mSUGRA scenario which 
consists of a point with an intermediate value of tan/? and a model line 
attached to it (SPS la) and of a "typical" mSUGRA point with relatively 
high tan/5 (SPS lb). The parameters for the SPS la point are ttiq = 100 
GeV, mi/2 = 250 GeV, Aq = -100 GeV, tan/5 = 10, /i > 0. For this point the 
Higgs masses are m^o = 113.7 GeV, m^o = 394.65 GeV, 112^0 = 394.9, mfj± 
= 403.6 GeV according to the Hdecay and Feynhiggsfast programs [20, 21]. 

The analysis is done for the center-of-mass energy of ^/s = 1 TeV, at 
which the cross section for the process e+e~ — *> HA is 2.5 fb. The luminosity 
assumed is 1000 fb~^. The branching ratio for the H (A) Higgs boson to bb 
is 0.64 (0.40), Ftoi = 0.785 GeV (1.251 GeV). 

The results presented in Fig. 16 are the mass sum and the mass difference 
after the selection cuts, kinematic fit and the final cuts for the mass difference 
and the mass sum respectively. The masses can be measured with precision 
of 1.3 GeV. The signal efficiency is 29% after selection cuts and 24% after 
cuts on di-jet mass sum and difference. The cross section can be measured 
with the relative uncertainty of 9%. 
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v/i [GeV 


(mH,mA) GeV] 


o-HA fb 




(150,100) 


33.62 




(200,100) 


25.30 




(250,100) 


16.61 


500 


(150,140) 


28.39 




(150,150) 


26.90 




(200,150) 


18.85 




(250,150) 


10.67 




(200,200) 


11.35 




(300,150) 


10.55 




(290,200) 


9.54 




(300,250) 


7.49 


800 


(300,300) 


5.70 




(350,350) 


2.23 




(400,150) 


6.46 




(400,200) 


5.17 




(400,250) 


3.70 



Table 1: Tree level cross sections o"ha for e'''e -^ HA expected for the Higgs 
boson mass hypotheses (r?7,H;''^A) considered in the study. Numbers are given 



for sin (/5 



a) 



-1. Cross sections are calculated including ISR at center-of- 



mass energies of -y/i = 500, 800 GeV. 
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Process 


ax BR fb] 


500 GeV 


800 GeV 


tt -^ w+bw-b 


6.69- 10^ 


1.65-102 


WW — > qqqq 


4.13 


10^ 


2.34 


10^ 


WW ^ qq£z/^ 


3.88 


10^ 


2.20 


10^ 


ZZ -^ qqqq 


3.11 


102 


1.74 


102 


ZZ ^ qq£+i- 


0.89 


102 


4.97 


101 


ZZ — > qqz/i^ 


1.78 


102 


9.96 


101 


qq(q = u,d, s,c,b) 


1.29 


10^ 


5.45 


103 


Wei^ -^ qqez/ 


5.08 


10^ 


8.40 


103 


Ze"'"e~ — »■ qqe+e~ 


4.46 


103 


4.10 


103 



Table 2: Topological cross sections a x BR of background processes at 
^ = 500, 800 GeV. 



Cut 


Tot. BG 


ZZ 


WW 


2-ferni. 


tt 


Eff.% 


Signal 




9380780 


157240 


2087383 


6799865 


336290 


100.0 


13618 


1 


6035870 


146909 


1957360 


3651766 


279839 


98.27 


13383 


2 


2675300 


105981 


1309760 


1047851 


211705 


91.34 


12438 


3 


1357080 


39789 


383774 


772600 


160915 


85.06 


11583 


4 


335995 


8167 


50257 


119448 


158123 


67.81 


9235 


5 


306623 


7832 


48523 


92288 


157981 


66.58 


9067 


6 


265119 


7499 


46254 


84886 


126480 


63.18 


8604 


7 


148085 


5975 


35270 


70218 


36621 


54.39 


7407 


8 


41912 


1761 


809 


10136 


29206 


53.22 


7247 


9 


1960 


272 


19 


473 


1196 


42.63 


5805 



Table 3: Number of remaining background events, signal efficiency and num- 
ber of expected signal events after each of the nine cuts for the mass point 
(mH,mA) = (150,100) in the HA -> bbbb channel at y/s = 500 GeV. The 
first column is the number of the cut, the second one is the number of total 
background events. 
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Cut 


Tot. BG 


ZZ 


WW 


2-ferni. 


tt 


Eff.% 


Signal 




4137520 


87211 


1183901 


2783826 


82579 


100.0 


3040 


1 


2376220 


69491 


922289 


1307850 


76592 


97.00 


2949 


2 


1007200 


41522 


462295 


432605 


70781 


93.98 


2857 


3 


542191 


15305 


154934 


319818 


52133 


85.67 


2604 


4 


90601 


2328 


21986 


42023 


24265 


82.22 


2499 


5 


78613 


2117 


20732 


31502 


24262 


80.29 


2441 


6 


54801 


1937 


16460 


29276 


7128 


71.31 


2168 


7 


12283 


534 


2075 


4430 


5244 


68.99 


2097 


8 


488 


50 


57 


175 


206 


44.92 


1366 



Table 4: Number of remaining background events, signal efficiency and num- 
ber of expected signal events after each of the eight cuts for the mass point 
(mH,mA) = (300,250) in the HA -^ bbbb channel at y/s = 800 GeV. The 
first column is the number of the cut, the second one is the number of total 
background events. 
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v^ [GeV] 


(mH,mA) [GeV 


Efficiency %] 


Number of events 


Tot. backgr. 


cuts 


2d. cut 


cuts 


2d.cut 


cuts 


2d. cut 




(150,100) 


43 


31 


5805 


4196 




132 




(200,100) 


41 


26 


4194 


2661 




129 




(250,100) 


36 


28 


2422 


1851 




182 


500 


(150,140) 


48 


39 


5469 


4518 


1960 


301 




(150,150) 


45 


41 


4944 


4407 




144 




(200,150) 


42 


27 


3192 


2036 




130 




(250,150) 


36 


22 


1534 


960 




156 




(200,200) 


37 


33 


1691 


1510 




185 




(300,150) 


43 


33 


1825 


1427 




38 




(290,200) 


46 


36 


1763 


1398 




53 




(300,250) 


45 


36 


1366 


1101 




13 


800 


(300,300) 


44 


35 


1011 


823 


488 


26 




(350,350) 


37 


31 


339 


277 




58 




(400,150) 


37 


26 


969 


678 




27 




(400,200) 


40 


28 


833 


584 




41 




(400,250) 


40 


28 


597 


414 




49 



Table 5: The signal efficiencies, the number of signal and total background 
events after selection cuts and after cuts on di-jet mass sum and difference in 
the HA -^ bbbb channel at ^/s = 500 and 800 GeV. The signal expectations 
are quoted for sin^(/3 — a) = 1 and the Higgs boson branching fractions of 
Br(H -> bb) = Br(A -^ bb) = 90%. 
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v^ [GeV] 


(mH,mA) [GeV] 


HA^l 


obr+r 


HA^- 


r+r^bb 


Tot.BG 


Eff. %] 


Events 


Eff. % 


Events 




(150,100) 


68 


1035 


73 


1110 






(200,100) 


67 


765 


72 


820 






(250,100) 


61 


460 


70 


525 




500 


(150,140) 


73 


940 


73 


940 


133490 




(150,150) 


74 


900 


74 


900 






(200,150) 


70 


595 


71 


605 






(250,150) 


63 


305 


70 


340 






(200,200) 


69 


355 


69 


355 






(300,150) 


60 


290 


72 


345 






(290,200) 


65 


280 


70 


305 






(300,250) 


64 


215 


67 


230 




800 


(300,300) 


64 


165 


64 


165 


52020 




(350,350) 


58 


55 


58 


55 






(400,150) 


52 


150 


71 


205 






(400,200) 


55 


130 


68 


160 






(400,250) 


56 


95 


66 


110 





Table 6: The HA -^ bbr+r" and HA 



T^T 



bb efficiencies, the number 



of signal and total background events after preselection at ^/s = 500 and 
800 GeV. The signal expectations are quoted for sin^(/3— a) = 1 and the Higgs 
boson branching fractions of Br(H — > bb) = Br(A -^ bb) = 90%, Br(H 
= Br(A -^ tt) = 10%. 



TT 
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v^ [GeV] 


(mH,mA) [GeV] 


HA^l 


obr+r 


HA^- 


r+r-bb 


Tot.BG 


Eff. %] 


Events 


Eff. % 


Events 




(150,100) 


42 


640 


44 


670 


5070 




(200,100) 


39 


445 


42 


480 


5890 




(250,100) 


21 


160 


28 


210 


2320 


500 


(150,140) 


41 


525 


42 


540 


3870 




(150,150) 


42 


510 


42 


510 


4275 




(200,150) 


34 


290 


38 


325 


2870 




(250,150) 


30 


145 


35 


170 


3815 




(200,200) 


29 


150 


29 


150 


2050 




(300,150) 


41 


195 


44 


210 


2023 




(290,200) 


41 


180 


43 


190 


1350 




(300,250) 


38 


130 


42 


140 


1056 


800 


(300,300) 


41 


105 


41 


105 


1609 




(350,350) 


40 


40 


40 


40 


2632 




(400,150) 


35 


100 


45 


130 


2528 




(400,200) 


37 


85 


45 


105 


2430 




(400,250) 


38 


65 


45 


75 


1842 



Table 7: The HA -^ bbr''"r~ and HA -^ r^r^bb signal efficiencies and the 
number of signal and total background events accepted into the final samples 
at a/s = 500, 800 GeV. The signal expectations are quoted for sin(/3 — a) = 
1 and the Higgs boson branching fractions of Br(H -^ bb) = Br(A -^ bb) = 
90%, Br(H -^ tt) = Br(A ^ tt) = 10%. 
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yi [GeV] 


(mH,mA) [GeV 


AniH [GeV 


AniA [GeV 


bbbb 


bbr+r 


Comb. 


bbbb 


bbr+r 


Comb. 




(150,100) 


0.12 


0.21 


0.10 


0.12 


0.18 


0.10 




(200,100) 


0.15 


0.45 


0.14 


0.15 


0.30 


0.13 




(250,100) 


0.27 


0.55 


0.24 


0.27 


0.46 


0.23 


500 


(150,140) 


0.13 


0.28 


0.12 


0.13 


0.25 


0.12 




(150,150) 


0.12 


0.40 


0.11 


0.12 


0.40 


0.11 




(200,150) 


0.26 


0.58 


0.24 


0.26 


0.40 


0.22 




(250,150) 


0.48 


1.40 


0.45 


0.48 


0.77 


0.41 




(200,200) 


0.31 


1.30 


0.30 


0.31 


1.30 


0.30 




(300,150) 


0.33 


0.44 


0.26 


0.33 


0.57 


0.29 




(290,200) 


0.33 


0.89 


0.31 


0.33 


0.76 


0.30 




(300,250) 


0.45 


1.08 


0.42 


0.45 


1.02 


0.41 


800 


(300,300) 


0.48 


1.34 


0.45 


0.48 


1.34 


0.45 




(350,350) 


0.97 


2.12 


0.88 


0.97 


1.12 


0.73 




(400,150) 


0.47 


1.36 


0.44 


0.47 


1.21 


0.44 




(400,200) 


0.74 


2.00 


0.69 


0.74 


1.88 


0.69 




(400,250) 


0.99 


2.70 


0.93 


0.99 


2.36 


0.91 



Table 8: Precision on the Higgs boson mass determination for the bbbb, 
bbr"'"r~ channels and their combination at y/s = 500 and 800 GeV. Num- 
bers are quoted for sin^(/5 — a) = 1 and Higgs boson branching fractions of 



Br(H -^ bb) = Br(A -^ bb) = 90%, Br(H -^ tt) = Br(A -^ tt) 



10%. 
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v^ [GeV] 


(mH,mA) [GeV 


Aa/ff [%■ 




HA -^ bbbb 


HA -^ bbr+r- 


HA -^ T+T-hh 




(150,100) 


1.6 


5.2 


4.9 




(200,100) 


2.0 


7.5 


5.8 




(250,100) 


2.4 


16.3 


8.9 


500 


(150,140) 


1.5 


7.2 


6.3 




(150,150) 


1.5 


4.2 


4.2 




(200,150) 


2.3 


9.7 


8.7 




(250,150) 


3.5 


22.6 


13.5 




(200,200) 


2.7 


8.1 


8.1 




(300,150) 


2.7 


10.0 


8.8 




(290,200) 


2.7 


11.2 


10.9 




(300,250) 


3.0 


13.8 


11.9 


800 


(300,300) 


3.5 


10.0 


10.0 




(350,350) 


6.6 


17.0 


17.0 




(400,150) 


3.9 


16.3 


11.6 




(400,200) 


4.3 


21.5 


19.8 




(400,250) 


5.2 


31.7 


28.5 



Table 9: Relative uncertainty in the topological cross section measurements 
for the bbbb, bbr"'"r~ channels at ^/s = 500 and 800 GeV. Numbers are 

and the Higgs boson branching fractions of 
= 90%, Br(H -^ tt) = Br(A -^ tt) = 10%. 



quoted for sin {(3 — a) = 
Br(H -^ bb) = Br(A -^ bb) 
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(mH,mA) GeV] 


rH,A [GcV 


STu [GeV 


ST A [GeV 


(150,100) 





<0.9 


<0.9 


5 


-U.6 
+1.3 


-U.6 
+1.3 


10 


-i.3 
+1.6 


-1.1 

+1.6 


(200,150) 





<1.7 


<1.8 


5 


-1.1 

+2.1 


-1.2 
+2.1 


10 


-2.U 
+2.6 


-2.U 
+3.0 


(300,250) 





<3.6 


<3.6 


5 


-2.5 
+2.9 


-2.5 
+2.9 


10 


-2.y 

+3.6 


-2.8 
+4.0 



Table 10: The accuracy of the Higgs boson width measurements in the case 
of large mass splitting between Higgs bosons. The upper limits on the Higgs 
boson width in the case of Fh = Fa = GeV correspond to 95% confidence 
level. 



(mH,mA) GeV] 


Th,a [GeV] 


5Tn,A [GeV] 


(150,150) 





<1.4 


5 


-l.U 

+1.4 


10 


-l.V 
+2.1 


(200,200) 





<3.3 


5 


-2.U 
+2.5 


10 


-3.2 
+3.7 


(300,300) 





<4.5 


5 


-3.2 

+3.7 


10 


-4.U 

+4.5 



Table 11: The accuracy on the Higgs boson width measurements in the case 
of Higgs boson mass degeneracy. The upper limits on the Higgs boson width 
in the case of Fh = Fa = GeV correspond to 95% confidence level. 



26 



^ 0.8 
0.6 








■ b-Net Brahms 
0.4 k c-N;et Brahms 

Q2 IQ brNetSimdet 
c-Net Simdet 



J I I I I I I I I I I I I I I I I I L 



0.2 0.4 0.6 0.8 1 

efficiency 



Figure 1: b-tag and the c-tag purity of the neural network (Section 4.1) versus 
efficiency curves for Z — > qq events simulated with SIMDET and BRAHMS 
at a center-of-mass energy of ^/s = 91.2 GeV. Filled squares (triangles) cor- 
respond to the neural network output for the b-events (c-events) simulated 
with BRAHMS. Open squares (triangles) correspond to the neural network 
output for the b-events (c-events) simulated with SIMDET. The neural net- 
works are trained on event samples simulated with SIMDET using the same 
variables and jet classification as in BRAHMS. 
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Figure 2: Distributions of the reconstructed invariant mass sum (upper figu- 
re) and mass difference (lower figure) in the sample of the HA — > bbbb events 
with {m-HymA) = (150,100) GeV at -y/i = 500 GeV. Dashed histograms show 
raw spectra obtained using only measured jet angles and energies. Solid 
histograms show spectra obtained after applying 4C kinematic fit. Three 
entries per event contribute to distributions, corresponding to three possible 
di-jet pairings. 
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Figure 3: The mass spectrum of the di-jet (upper figure) and di-tau (lower 
figure) systems in the sample of the HA -^ bbr"'"r~ events with (mH,''^A) = 
(150,100) GeV at a/s = 500 GeV. Dashed histograms represent distributions 
before applying kinematic fit, solid - after applying kinematic fit. 
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Figure 4: Distributions of the selection variables (visible energy, number of 
tracks per jet A^fracfcs/jet, cos6t and a thrust value) in the HA — > bbbb chan- 
nel with (mH,mA) = (150,100) GeV at ^/s = 500 GeV. The signal distribu- 
tions are shown with arbitrary normalisation. The vertical lines and arrows 
indicate cuts, imposed on these variables. The distributions are shown after 
all cuts preceding the current variable. 
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Figure 5: Distributions of the selection variables (jet resolution parameters 
logio ?/34 ^^^ logio 1/56 ^T^^ b-tag Variables -B12 and -B34) in the HA — ^ bbbb 
channel with (17111,171^) = (150,100) GeV at ^/s = 500 GeV. The signal distri- 
butions are shown with arbitrary normalisation. The vertical lines and arrows 
indicate cuts, imposed on these variables. The distributions are shown after 
all cuts preceding the current variable. 
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Figure 6: Variables used to construct signal likelihood (thrust value, cos 6t, 
visible energy E^g and number of energy flow objects Nenfiow) in the 
HA -^ bbr+r^, r+r^bb channels with {mi{,mA_) = (150,100) GeV at a/s = 
500 GeV. The distributions are shown after cuts 1-6 (Section 5.2). 
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Figure 7: Variables used to construct signal likelihood (minimal opening 
angle between any two jets ajetjet and b-tag variables Bi and B2) in the 
HA -^ bbr+r^, r+r^bb channels with (mH,"7,A) = (150,100) GeV at ^/s = 
500 GeV. The distributions are shown after cuts 1-6 (Section 5.2). 
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Figure 8: Distribution of Lha in the HA -^ bbr+r^, r+r~bb channels with 
(mH,mA) = (150,100) GeV at ^i = 500 GeV. The vertical line and arrow 
indicate cut placed on the likelihood. 
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Figure 9: Distributions of the di-jet mass sum (upper figure) and di-jet mass 
difference (lower figure) in the HA — ^ bbbb channel for Higgs boson mass 
hypothesis (mH,rraA) = (300,250) GeV at y/s = 800 GeV after selection cuts 
and kinematic fit. The three components of the signal are shown separately: 
two combinatorial parts and one of the real signal. The background from the 
other processes is presented as well. 
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Figure 10: Upper figure: distribution of tlie di-jet mass sum after selection 
cuts, kinematic fit and cut on di-jet mass difference. Lower figure: distribu- 
tion of tlie di-jet mass difference after selection cuts, kinematic fit and cut on 
di-jet mass sum. Both distributions are in the HA -^ bbbb channel for the 
Higgs boson mass hypothesis {mii,m\) = (300,250) GeV at y/s = 800 GeV. 
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Figure 11: Distributions of di-jet mass and di-tau mass in the 
HA -^ hhr^T^ ,T^T~hh channels for the case of {mii,mA) = (200,150) GeV 
at v^ = 500 GeV. 
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Figure 12: Distributions of di-jet mass and di-tau mass in the 
HA -^ bbr+r^, r^~r~bb channels for the case of {mE_,mA) = (300,300) GeV 
at ^ = 800 GeV. 
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Figure 13: Discovery significance as a function of the Higgs boson mass 

(assuming rriH = m^i) in the HA —>■ bbbb channel at ^/s = 800 GeV (upper 

figure 

(lower figure) where t]'^ is assumed e"'"e 



7] as a function of the Higgs boson mass for the 5a discovery limit 



for sin {/3 — a) 
assumed to be 90%. 



* HA cross section relative to that 
1. Branching ratios Br(H -^ bb) and Br (A -^ bb) are 
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Figure 14: The reconstructed di-jet mass spectrum in the HA — »• bbr"''r~ 
sample (upper figure) and the reconstructed di-tau mass spectrum in the 
HA -^ r+r~bb sample (lower figure) for the different Higgs boson widths Fh 
= 0, 5, 10 GeV for m^.A = 200 GeV at ^i = 500 GeV. 
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Figure 15: The log-likelihood as a function of Fh and Fa for the case of 
Th = Fa = 5 GeV and (mH,mA) = (200,150) GeV at ^i = 500 GeV. 
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Figure 16: Upper figure: distribution of tlie di-jet mass sum after selection 
cuts, kinematic fit and cut on di-jet mass difference. Lower figure: distribu- 
tion of the di-jet mass difference after selection cuts, kinematic fit and cut 
on di-jet mass sum. Both distributions are in the HA -^ bbbb channel for 
the SPS la benchmark point with Higgs boson mass hypothesis (mH;''^A) = 
(394.90,394.65) GeV at v^ = 1 TeV. 
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